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H ■ Abstract 

A multi-scale characterization of the field concentrations inside composite and poly- 
crystalline media is developed. We focus on gradient fields associated with the intensive 
quantities given by the temperature and the electric potential. In the linear regime 
these quantities are modeled by the solution of a second order elliptic partial dif- 
ferential equation with oscillatory coefficients. The characteristic length scale of the 
heterogeneity relative to the sample size is denoted by e and the intensive quantity is 
| denoted by u £ . Field concentrations are measured using the LP norm of the gradient 

field ||Vn e || i p (D) for 2 < p < oo. The analysis focuses on the case when < £ < 1. 
Explicit lower bounds on liminf e ^o || Vu e \\ L v( D) are developed. These bounds provide a 
way to rigorously assess field concentrations generated by the microgeometry without 
having to compute the actual field u e . 
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g ! 1 Introduction 

The initiation of failure inside heterogeneous media is a multi-scale phenomena. Loads ap- 
plied at the structural scale are often amplified by the microstructure creating local zones 
of high field concentration. The local amplification of the applied field creates conditions 
that are favorable for failure initiation ^U|. This paper focuses on gradient fields associated 
with the intensive quantities given by the temperature and the electric potential inside het- 
erogeneous media. The local integrability of the gradient directly correlates with singularity 
strength which influences the onset of failure such as dielectric breakdown. 

In this work it is shown how to assess the L p integrability of the gradient fields in mi- 
crostructured media by investigating the multi-scale integrability of suitably defined quan- 
tities. The analysis is carried out with minimal regularity assumptions on the coefficients 
describing the local properties inside the heterogeneous media. The results are described in 
terms of the p th order moments of the solution of two-scale corrector problems. The quan- 
tities are sensitive to microscopic field concentrations and can become divergent for p > 2. 
This is in contrast to the well known effective constitutive properties which are based upon 
local averages and are bounded above independently of the microgeometry. 
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The results given here are presented in the context of two-scale homogenization pQ, |2Uj . 
We consider a bounded domain Q in R n , n > 2. A common microstructure that admits a 
two-scale description is a simple generalization of a uniformly periodic microstructure and is 
described as follows. Consider a partition of the domain Q made up of measurable subsets Qg, 
£ = 1,2, . . . , K such that Q = uf =l VLg. Inside each subdomain Qg we place a different periodic 
microstructure made from N anisotropic heat conductors. This type of microstructure will 
be referred to as a piece wise periodic microstructure. Well known engineering composites 
that are modeled by piecewise periodic microstructures include multi-ply fiber reinforced 
laminates 0, [21] and [23! • 

The thermal conductivity tensor for the piecewise periodic microstructure is described as 
follows. The indicator function for each of the subdomains fig is denoted by Xn,( x )> taking 
the value 1 for points in Qg and zero outside. In order to describe the periodic microstructure 
inside the I th subdomain we introduce the unit period cell Q. The configuration of the N 
phases inside Q is described by the indicator functions x\(y)i i — 1, ■ ■ ■ ,N associated with 
each phase. Here Xg{y) = 1 f° r points inside the i th phase and zero outside. The length scale 
of the microstructure relative to the size of the domain Q is given by e k — 1/k, k = 1, 2 . . .. 
The microstructure is obtained by rescaling the configuration inside the unit period cell. The 
indicator function of the i th conductor in the microstructured composite is given by 

K 

X- fc (x) = Xi(x,x/e fc ) = ^Xo,(x)xK x M)- (1.1) 

i 

The local conductivity tensor A £k has a two-scale structure and is given by 

TV 

A £ *(x) = A(x,x/e k ) = ^A l Xl (x,x/e fc ). (1.2) 



Other heterogeneous media that are amenable to similar or more general two-scale de- 
scriptions include polycrystalline materials such as metals and ceramics. We state the general 
hypotheses under which the two-scale homogenization theory applies, see and [2]. It is 
assumed that A(x, y) is a matrix defined on Q x Q and there exist positive numbers a < (3 
such that for every vector 77 in R 3 that 

a\ V \ 2 <A(x,y)<(3\ V \ 2 . (1.3) 

The conductivity Aij(x, y) is Q-periodic in the second variable, such that Aj( x , x/e fc ) is 
measurable, satisfies 



lim 

e fc ->0 



n 



x 2 
Aij(x., — ) 



dx = I |Ay(x,y)| 2 dxdy (1.4) 



IflxQ 

and for any suitable two scale trial field ip(x,y) that 



lim / Ay(x, — )ip(x,x/e k )dx = / A ij (x,y)i/;(x,y) dxdy. (1.5) 
£ k^°Jn £k Jnxo 
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The convergence given by ()1.5|) is a weak convergence and is known as two-scale convergence 
[T] (20]. The space of suitable two-scale trials is denoted by L 2 [D; C per (Q)]. Here C per {Q) 
denotes Q-periodic continuous functions defined on R 3 and the space L 2 [D; C per (Q)] is the 
space of functions h : Q — > C per (Q) which are measurable and satisfy J n ||/i(x)||^ er(Q) c?x < oo. 
The norm ||/i(x) ||c per (Q) ^ s defined by sup yg g |/t(x,y)|. In what follows no other regularity 
hypothesis on the conductivity matrix A(x, y) is made. 

The temperature field u £k associated with the conductivity tensor field A £k (x) = A(x, x/e^) 
is the solution of the equilibrium equation 

- div (A £k (x)Vu £k ) = / in Q (1.6) 

with the boundary conditions given by u 6k = on dflo and n ■ A £ky \7u 6k = g on 8Qn with 
dtt = dVL D U dVL N . 

In what follows we consider the limit as tends to zero. We fix a subdomain D of Q 
and derive lower bounds on 

liminf ||VM £fc || L P (D) . (1.7) 

The lower bound is expressed in terms of a two-scale integral that encodes the field 
amplification properties of the microstructure. It is formulated in terms of the solution of 
the homogenized problem together with a local corrector matrix that captures the interaction 
between the periodic microstructure and the gradients of the homogenized temperature field. 
The bounds introduced here provide a rigorous way to assess field concentrations generated 
by the microgeometry without having to compute the full solution u £k . 

The lower bound is given in terms of the solutions w/(x, y) to the local periodic problem. 
For each x in Q the function u/(x, y) is a Q periodic function of the second variable y and 
is a solution of 

div y (A(x,y)(VX(x,y) + e i )) =0, (1.8) 
with f Q -u/(x, y) dy = 0. The corrector matrix P(x, y) is defined by 

P tJ (x,y) = 9 2/j ^(x,y) + e}. (1.9) 
The associated effective conductivity tensor A E (x) is given by 

A E (x)= [ A(x,y)P(x,y)dy. (1.10) 
J Q 

The two-scale homogenization theory gives the following theorem pQ. 
Theorem 1.1. Two-scale Homogenization Theorem 

The sequence of solutions {u £k } £k> o of M.b]) converges weakly to u H (x) in P 1 (fi) where u H 
is the solution of the homogenized problem 

- div (i B (x)V/(x)) = /(x), inil, 

u H (x) = 0, on dn^, and 
n-A E Vu H = g, ondn N . (1.11) 
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The field concentration functions of order p are defined by 

i/p 



f p (x,Vu H (x)) = Qf |P(x,y)Vn H (x)rdy^ P , 2 < p < oo (1.12) 

and / p (x, V« fl (x)) < / g (x, Vm f (x)) for p < q. It is clear that / p corresponds to a p t/l order 
moment of the corrector matrix (jl.9|) and 

/oo(x,Vu ff (x)) = limf/ |P(x,y)V^(x)| p C /yV P . (1.13) 

Theorem 1.2. Lower Bounds on Field Concentrations 
For 2 < p < oo 

(7 (/ P (x,V^(x))) p dx^ /P < liminf ||Vn £fc || LP(D) . (1.14) 

For multi-phase conductivity problems with coefficients described by (jl.2|) the field con- 
centration functions of order p are defined for each phase and are given by 

/;(x,V« H (x))= Xi(x,y)|P(x ) y)V« H (x)p t iy^ " , z = l,...,JV, 2 < p < oo (1.15) 

and /p(x, Vw fl (x)) < /*(x, Vw fl (x)) for p < q. As before one defines 

/^(x,W(x)) = lim (/ Xj (x,y)|P(x,y)V^(x)rrf y y P . (1.16) 



For this case lower bounds on 



liminf \\x £ i k Vu £k \\ LP{D) (1.17) 



are given by the following theorem. 

Theorem 1.3. Lower Bounds for Multi-phase Composites 
For 2 < p < oo 

i/p 



Qf (/;(x,V^(x))) P dx^) ' < liminf|| X ^Vn £fc || LP(D) . (1.18) 



The bounds can be applied to develop a Chebyshev Inequality for the distribution func- 
tions associated with the sequence {Xi fc |VM £fc |} £fc >o- Here the distribution function A^ fe (.D,i) 
gives the measure of the set inside D where ^*|Vit efc | > t. 

Arguing as in Proposition 2.1 of J3| and combining with ()1.18|) gives the following 
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Theorem 1.4. Homogenized Chebyshev Inequality 



limsupAf^At) <r p (/ fl (/;(x,Vm h (x))) P (Ix) < rHiminf Vw 




LP(D)- 



(1.19) 



It is pointed out that Theorems 11.21 and 11.31 are obtained using the minimum regularity 
assumptions on the coefficients A £k . Because of this the hypotheses of (Theorem 2.6, pQ) do 
not apply and one can not take advantage of the strong convergence given in that Theorem. 
Instead the theorems are proved using a perturbation approach introduced in see 
Section 2. 

The lower bounds are sensitive to the presence of singularities generated by the mi- 
crostructure. To illustrate this we consider a microstructure made from a periodic distribu- 
tion of uniaxial crystallites embedded in an isotropic matrix of unit conductivity. The period 
cell for the composite is illustrated in Figure 1. Each crystallite occupies a sphere and has 
conductivity Ai in the radial direction and A2 in the tangential direction. The dispersion of 
the N crystallites is specified by Uf B(y l , 77) where B(y e , rg) denotes the £ th sphere centered 
at y e with radius rg. Each crystallite has a conductivity tensor given by 



where n = (y — y £ )/|y — y e \ for y in I?(y £ ,r^) and I is the 3x3 identity. Outside the 
crystallites we set A(y) = I. It is supposed that the aggregate of crystallites occupy an 
area fraction < 9 < 1 of the unit period cell. It is noted that the conductivity inside each 
crystallite is precisely the one employed in the Schulgasser sphere assemblage [2*4] . 

When a constant gradient field is applied to a single isolated crystallite and Ai > A 2 the 
crystallite exhibits a gradient field singularity at its center. In what follows we use the lower 
bound ()1.14|) to show how this local information effects the integrability of the sequence 
{Vw efe } £fe >o- We form A 6k = A(x/ek) and consider solutions u 6k of Ijl.fij) . To fix ideas we 
choose / to be in L r (Q) for r > 3 and g to be in L 2 {dVt^). In what follows A2 is restricted to 
lie in the interval 1/2 < A2 < 1 and Ai = 1/ (2A2 — 1). For this choice it is shown in Section 
3 that the homogenized temperature field u H is the solution of (|1.11|) with A E = I. 

For D compactly contained in Q it follows from the LP theory ^7j, that || Vu h \\lp(d) < 00 
for every 1 < p < 00. On the other hand calculation and application of Theorem 11.21 shows 
that 



A(y) = Ain <g> n + X 2 (I - n ® n) 



(1.20) 



LB{p) x ||ViT 



< liminf || Vtt' 



(1.21) 



where 




for, p < 
for, p > 



3 



3 



(1.22) 
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For a fixed choice of A2 the value p c = 2 (i-\ 2 ) satisfies 3 < p c < +00 and 

liminf || \7u ek \\ LP(D) = +00, for p > p c . (1-23) 

This is in stark contrast to the L p integrability of the the gradient of the homogenized solution 
which holds for any p < +00. It is clear for this example that the information carried by 
the homogenized problem is not adequate and misses the singular behavior exhibited by the 
sequence {Vu £k } £k> Q. This example shows that failure initiation criteria based solely upon 
the solution of the homogenized equations will be optimistic. The inequalities given above 
are established in Section 3. 

The maximum integrability exponent for the gradient of the solution of the local problem 
(jl.Sj) is referred to as the threshold exponent for the composite. The threshold exponent is 
introduced in the work of Milton ^H] and measures the worst singularity of the gradient field. 
The threshold exponent is precisely p c for the local problem considered here and corresponds 
to the divergence in the lower bound for p > p c . 




Figure 1: Unit period cell with Schulgasser crystallites embedded inside a material with unit 
thermal conductivity. 

Next we consider an example for which the sequence {Vw £fc } efe >o is uniformly bounded 
in LP for some class of coefficients and right hand sides /. For this case we show that the 
lower bound given in Theorem 11.21 is attained. In this example we make use of the a priori 
estimates for {\7u £k } £k> o developed in the Theorem 4 of Avellaneda and Lin Let Q be 
a C l ' a domain (0 < a < 1) and suppose for < 7 < 1, C > 0, that A(y) G C 7 (R n ) and 
||v4(y) ||c7(r«) < C. Then we choose A 6k = A(x/£&). For 5 > suppose 2 < q < n + 5 and 
/ G L q and set 1/q = l/q—l/(n + 8). Given these choices we consider the Wq 1 ' 2 ^) solutions 
u £k of 

- div (A £k (x.)Vu £k ) = f in ft. (1.24) 

It is shown in Section 4 that ()1.14|) holds with equality for every p such that p < q. In fact 
it is seen more generally that for p < q and any Caratheodory function ip : D x R n — > R 
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satisfying 

|^(x,r?)| < \v\ P , for a.e. x G D and r] G R 3 , (1.25) 

that 

lim / ^(x, Vw £fc ( x )) rfx = / / ^(x,P(y)VM^(x))rfyrfx. (1.26) 
£ k^J D JdJq 

This is established in Section 4. 

It is anticipated that there are several classes of conductivity coefficients and right hand 
sides / for which the lower bounds are attained. In this direction we point out the the recent 
higher regularity results given in 0, [0], [TT] . [TJ] and [23] . 

We conclude noting that the analogues of the field concentration functions ()1.12|) and 
()1.15|) have appeared earlier in the contexts of G-convergence and random media, see 
and [Hj. In those treatments they are shown to provide upper bounds for the distribution 
function of the local stress and electric field for G-convergent sequences of elasticity tensors 
and random dielectric tensors. 



2 Derivation of the lower bounds 

We recall the weak formulation of the Ek > problem given by (jl.fij) . Let V denote the 
closure in if x (f2) of all smooth functions that vanish on dflo- We suppose that / is in L 2 (Q) 
and g belongs to L 2 {dVL^). The function u £k belonging to V is the solution of the weak 
formulation of the boundary value problem given by 

/ A(pc,x/e k )Vu Sk ■ Vy?cix = / ftpdy.+ I gipds, (2.1) 
Jn Jn Jan N 

for every ip in V. Here ds is an element of surface area. 

In order to express the two-scale weak formulation of we introduce the following 

function spaces. The space of square integrable Q-periodic mean zero functions with square 
integrable derivatives is denoted by Hp er (Q)/H. The norm of an element v in this space 
is denoted by IHIh^cqvr- The space of measurable functions h from Q to Hp er (Q) /R 
for which f n ||/i(x)||^i /R dx < oo is denoted by L 2 [Q; Hp er (Q)/H}. This function space 

was introduced for the description of the two-scale homogenized problem in j201 • The weak 
formulation of the two-scale homogenized problem (jl.llj) is given by the unfolded variational 
principle PQ, [7], HE]. 

Theorem 2.1. Unfolded Variational Principle 

The pair (u H ,ui) is the unique solution in V x L 2 [Q; Hl er (Q)/'R] of 

[ [ A(x,y)(V/(x) + V y « 1 (x,y))- (V^(x) + V y <^i(x, y)) rfy rfx 
Jn Jq 

f(pdx+ / gpds, (2.2) 
< > Jdn N 
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for every (<p,<pi) in V x L 2 [Q; Hp er (Q)/H]. Moreover 

Vu H + V yMl (x, y) = P(x, y)Vu H (x). (2.3) 



In order to establish Theorems 11.21 and 11.31 we recall the function spaces used in the 
description of two-scale convergence [TH]. The space C per (Q) denotes Q-periodic continu- 
ous functions defined on R 3 . For 1 < r < oo the space L r [D; C per (Q)] is the space of 
functions h : D — > C per (Q) which are measurable and satisfy f D ||/i(x)||£, (Q) dx < oo. 
Here \\h(x) \\c per (Q) = su P y eQ IM X > y)|- The intersection of the spaces L°°(D x Q) and 
L r [D]C per (Q)} is denoted by by V r . For 1 < r < oo we introduce 1 < r' < oo such 
that - + ^7 = 1. We establish Theorems 11.21 and II .31 with the aid of the following Lemmas. 

Lemma 2.1. Localization Lemma 

Fix a domain of interest D inside Q. Let q(x, y) be any test function in V r then one can 
pass to the limit e & — > in the sequence of solutions {u Sk } ek>0 of M.b}) to obtain: 

lim / g(x,x/ £fc )|Wf dx = I [ g(x,y)|P(x,y)V^(x)| 2 rfyrfx. (2.4) 
£ k->oj D J D Jq 



For multi-phase composites with coefficients described by (jl.2j) we restrict attention inside 
each phase and state the following lemma. 

Lemma 2.2. Localization Lemma in Multi-phase Composites 

Let g(x, y) be any test function in V r then one can pass to the limit — > in the sequence 
of solutions {u £k } Ek> o of M.b}) to obtain: 

lim / g(x, x/e fe ) X; fc (x) \Vu £k \ 2 dx 

= 11 g(x,y) X ,(x,y)|P(x,y)V^(x)| 2 dydx. (2.5) 
J d Jq 



The proofs of Lemmas 12.11 and 12.21 are given at the end of this section. 

To illustrate the ideas we use Lemma f2. 21 to establish Theorem II .31 noting that Theorem 
II. 21 follows from Lemma 12. II in the same way. 

Proof of Theorem 11.31 For each > we apply Holder's inequality to the left side 
of to obtain 

/ / g(x,y)x,(x,y)|P(x,y)Vn H (x)| 2 rfyrfx 
Jd Jq 

<\im(J | 9 (x,x/e fc )| r da^ ' lim inf (J x?(x) I Vw £fe | 2r ' dxj ' . (2.6) 
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Noting PH that 

i™ (/, k( x > x / £ fe)r^ = (/d Jq k( x ; y)l rrf y rfx ) 7 = lk( x ) y)IU'-(DxQ) (2.7) 



we obtain 



Id Iq 9( x >y)Xi( x >y) |P(x,y)Vu H (x)| 2 dydx 
||g(x,y)|| £ r (Dx0) 



< liminf^y Xi fc (x)|VM £fc | 2r 'c/xV(2.8 



Since is dense in Z/(.D x Q) we take the supremum of the left hand side of ()2.8j) over V r 
to find that 

jf y Xi (x,y) |P(x,y)V^(x)| 2 ^ydx) 7 < limmf Qf x?*(x)|Vu e M V <*x) 7 (2.9) 

Theorem 1 1 . HI follows for 2 < p < oo upon taking the square root on both sides of (|2.9j) . The 
case p = 2 follows immediately upon choosing g(x, y) = 1 in Lemma [2.21 

We conclude by providing the proof of Lemma 12.21 and note that the proof of Lemma 12.11 
is identical. 

Proof of Lemma 12.21 The indicator function of the set of interest D is denoted by 
Xd(x). We choose a test function g(x, y) in V r and set p(x, y) = XD( x )Xi( x > y)<?( x ? y)- F° r 
5(3 > we form the perturbed conductivity tensor Ay(x, y) = Ajj(x, y) + 5/3p(x, y)<%. We 
choose 5/5 sufficiently small so that A(x, y) satisfies (jl.3j) . By construction A(x, x/e^) is 
measurable and satisfies (jl.4j) and (jl.5j) . Consider the associated solution £t £fc in V of the 
weak formulation of the boundary value problem given by 

/ A(x, x/£fc)V£t £fe • V<£> <ix = / f(pdx+ / g^fds, for every <£> in V. (2.10) 

Set u £k = u £k + <5w £fc and subtraction of ()2.1|) from ()2.10|) gives 

/ i(x,x/e fe )V5u £fc ■ V^rfx+ / 6Pp(x,x./e k )Vu £k ■ V^dx = 0. (2.11) 

Choosing <^ = u 6k in ()2.11|) and application of the identity 

A(x,x/e k )Vu £k ■ V5u £k dx = [ fSu £k dx+[ g5u £k ds, (2.12) 
<» in Jdn N 



gives 



5(3 x i p{x,x/e k )\Vu £k \ 2 dx + T £k = - f5u £k dx- g5u £k ds, (2.13) 
Jn Jn Jdn N 



where 



5(3 x f p(x,x/e k )(V5u £k ) ■ Vu £k dx. (2.14) 
Jn 



Next set ip = 5u £k in ()2.11|) and it follows from Cauchy's inequality and ()1.3|) that 

\\V8u £k \\ L 2 (n) <C6(3 } (2.15) 

where here and throughout C denotes a generic a constant independent of e&. From this it 
is evident that 

\T £k \<C5(3 2 . (2.16) 

Next we pass to the — > limit and apply Theorems 11.11 and 12.11 to find that the 
sequence {u £k } £k>0 converges weakly in to u H , where (u H ,ui) is the solution in 

VxL 2 [fi;^ er (Q)/R]of 

I(x, y)(Vu H (x) + V y «i(x, y)) ■ (V^(x) + V y ^(x, y)) dy dx 




'n 



/y>dx + / gtpds, (2.17) 

for every (<p, <pi) in V x L 2 [f2; ifp er (Q)/R]. Set u H — u H = 5u H , U\ — U\ = 8u\ and subtraction 
of (JUJ) from lj2T7|» gives 




'n 



+ 



A(x,y)(V«J« ff (x) + V y <J«i(x,y)) • (Vy?(x) + V y ^(x, y)) dy rfx 

/ / 5/3p(x,y)(V^(x)+V y «i(x,y))-(V^(x)+V y ^(x,y))dydx = 0.(2.18) 
</q Jq 

Choosing (<p,ipi) = {u H ,u\) in ()2.18|) together with the identity 

/ / A(x,y)(VM H (x) + V y M 1 (x,y))-(V^(x) + V y 5M 1 (x,y))c/yc/x 
Jn Jq 

f5u H dx+ [ g5u H ds. (2.19) 



gives 




5(3 x p(x,y)|P(x,y)V^(x)| i rfyrfx + T 



= - / fSu dx - / g5u H ds, (2.20) 

Jn JdQ N 

where 

f = 8/3x f f p(x,y)(V(5^ + V y (5 Ml (x,y))-(V^ + V yM i(x,y))rfx. (2.21) 
Jn Jq 

Next set (<p, <pi) = (5u H , in ()2.18|) and it follows from Cauchy's inequality and ()1.3|) that 

|| + V^xU^xq) < (2.22) 
10 



and it follows easily that 

\f\ < C5f3 2 . (2.23) 

Taking the — > limit in ()2.13|) noting that \im ek _ 5u £k = 5u H (weakly in H 1 ^)) and 
recalling (j2.16|) gives 

6px\im I p(x,x/£ fc )|VM £fe | 2 c/x + 0(5/? 2 ) = - I f5u H dx- [ g5u H ds. (2.24) 
£ fc^° Jn Jan N 

Lemma l2~2l now follows immediately from ()2.20)1 . ()2.23|1 and ()2.24|1 and identifying like powers 
of 5(3. 

3 Explicit lower bounds for aggregates of Schulgasser 
crystallites 

In this section we derive the lower bound (jl.22|) for the microstructure consisting of Schul- 
gasser crystallites embedded within a homogeneous matrix with unit thermal conductivity. 
The temperature field inside the unit period cell $*(y) = w l (y) + yi is the solution of the 
local problem 

div y (A(y)(V y «/(y) + e l ))=0, (3.1) 

with w l Q-periodic and Jq w l (y) dy = 0. For this microstructure A(y) is given by (jl.2U|) for 
y in B(y e ,ri) and A(y) = I outside. Here we restrict A 2 to the interval (1/2, 1) and choose 
Ai so that Ai = 1/(2A 2 — 1). A calculation shows that the solution $ ? (y) is given by 

*<_/ y<. yeQ\uf =1 5(y £ ,r,), 

' ~ I r}- a \y - y^- 1 ^ - y|) + yf, y G 5(yV £ ) ' ^ 

where a = 2X 2 — 1. The corrector matrix P(y) is given by 

P(v) _f I, yeQ\uf =1 £(yV,), 

[y) I r]- a |y-y^- 1 (/ + (o ! -l)n®n) ) yeB(y / l r / ) ' [6 ' 6) 

where n = (y — y^)/|y — y e \ for y G B(y e ,rg). A direct calculation shows that 

A E = [ A(y)P(y)dy = I. (3.4) 
J Q 

Next we provide the lower bound for J n |P(y)Vw^(x)| p <iy dx. Note for any rj in R 3 that 
P T (y)P(y)r]-r] = \P{y)ri\ 2 and the smallest eigenvalue A(y) of P T (y)P(y) delivers the lower 
bound A(y)|f?| 2 < \P(y)r]\ 2 and 

/ / \(y) p/2 \Vu H (x)\ p dydx< [ [ \P(y)X7u H (x)\ p dydx. (3.5) 
Jn J q Jn Jq 

Calculation shows that 

A(y)=a 2 rf- Q) |y-y? (Q - 1) (3-6) 

for y G B(y e ,ri) and A(y) = 1 for y G Q \ UfB(y e , rg). The lower bound (jl.22|) follows 
upon substitution of ()3.6|) into (j3.5jl . 
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4 Optimality of the lower bounds 

Conditions are presented on / and A(y) for which the lower bound ()1.14j) is attained for 
a range of exponents 2 < p < q. We suppose as in Avellaneda and Lin |H] that Q is a 
C 1 ' domain (0 < a < 1) and suppose for < 7 < 1, < C, that A(y) e C 7 (R n ) and 
||A(y)||c7(Rn) < C. We set A £k = A(x/s k ). For 5 > suppose 2 < g < n + 5 and f £ L q 
and set 1/g = l/q — l/(n + 5). Given these choices we consider the Wq 1 ' 2 ^) solutions u £k of 

- div (A £k (x.)Vu £k ) = f in n. (4.1) 

Theorem 4 of [3] shows that there exists a constant independent of e k for which 

||V^|| £4(n) < C||/|| w(n) (4.2) 

holds for every Sk > 0. Subject to these hypotheses it will be shown that the lower bound 
()1.14j) is attained for p < q. 

Passing to a subsequence if necessary we start by considering the Young measure v 
associated with the sequence {P(x./ek)'Vu H (x)} ek>0 . Here v is represented by a family 
of probability measures v = {z/ x } xg Q depending measurably on x. We denote the set of 
continuous functions (p defined on R" such that lim^oo <p(i]) = by Co(R n ). Elementary 
arguments show that 

<v x ,(p>= / (p(r))dv x (r]) = / (p(P(z)\7u H (x.))dz, a.e. x£!l (4.3) 
Jit" Jq 

for every ip in Co(R n ). From corrector theory (T^j there exists an exponent r > 1 for which 
one has the strong convergence 

lim \\Vu £k - P(x/e k )Vu H \\ L r {n) = 0. (4.4) 

The strong convergence ()4.4j) shows that both sequences {Vu £k } ek> o and {P(x/e / t)Vn i/ (x)} efe>0 
share the same Young measure see for example Lemma 6.3 of [221 ■ From ()4.2|) it follows on 
passage to a subsequence if necessary that {|VM £fc | p } £fc is weakly convergent in L 1 (f2) thus 

lim f \Vu £k \ p dx= [ [ \rj\ p du x (rj) dx = [ [ \P(z)Vu H \x)\ p dzdx, (4.5) 

£k^0J D J D J D Jq 

and optimality follows. Last, it follows immediately from Proposition 6.5 of (221 that for 
every Caratheodory function ip('x,r)) satisfying the growth condition (jl.25|) that (on passage 
to a subsequence if necessary) 

lini / ^(x, Vu £k ) dx = / / ^(x, rj) dvjjj) dx (4.6) 



D JD JH 



and (jl.26|) follows since ()4.3|) implies that 



ip(x,r])dv x (T})d]c= / / ijj(x,P(z)Vu H (x.))dzdx. (4.7) 

D JH n JD JQ 
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